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Employing the maximum entropy method we extract the spectral functions from meson correlators at four 
lattice spacings in quenched QCD with the Wilson quark action. We confirm that the masses and decay constants, 
obtained from the position and the area of peaks, agree well with the results from the conventional exponential 
fit. For the first excited state, we obtain m^^ = 660(590) MeV, nip^ = 1540(570) MeV, and fp^ = 0.085(36) in 
the continuum Hmit. 



1. Introduction 

The spectral function f{uj) of hadrons includes 
information such as masses and decay constants 
for various bound states and the continuum spec- 
trum for multi-particle states. Recently the max- 
imum entropy method |Q| (MEM) has been em- 
ployed in order to extract the spectral function 
from lattice QCD data ||,|. 

In this article we apply the MEM to the 
quenched lattice QCD data for pseudoscalar(PS) 
and vector (V) mesons previously calculated at 
/? = 5.90, 6.10, 6.25, and 6.47 g. We present 
the corresponding spectral functions for the case 
of a point source, and check the reliability of the 
MEM by comparing the masses and decay con- 
stants from the spectral function with those from 
the exponential fit. We then extract the masses 
and decay constants for the first excited state 
in the continuum limit. For details we refer to 
ref. §. 

2. MELximum entropy method 

In the MEM the spectral function is determined 
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by maximizing the quantity P[F\DH], which is 
the conditional probability of the spectral func- 
tion F for a given data D and all prior knowl- 
edge H such as f{u!) > 0. By Bayes's theorem 
this probabihty is replaced with F[D\FH]F[F\H], 
where P[D|Fi/] is proportional to exp(— x^/2) 
with the standard chi-squared, and P[i^|iJ] 
to exp(Q!S'(/)) with a a real positive parameter. 
Here the entropy S{f) ||l| is defined by /(w) and a 
positive function called model. Combining these 
two factors, wc can determine the most probable 
spectral function as the solution of the equation 
dQa/df = 0, where Q„ = aS{f) - xV2. Finally 
the result is averaged over a with a weight factor 
P[a\DH]. 

3. Results 

In Fig. H a typical result for the spectral func- 
tions of p meson is presented at /3 = 5.90 and 
6.47 for three different values of mT^/mp. The 
parameters of the MEM analysis are compiled in 
ref. |]. At each p one observes, as expected, that 
the peaks move to smaller Lua as quark mass de- 
creases. As /3 increases, the peaks of the ground 
and the first excited state also move to smaller 
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Figure 1. Spectral functions at /3 = 5.90 and 6.47 
for p meson. 

uja as expected. In addition the number of peaks 
increases, since higher excited states appear be- 
low the cutoff 7r/a at larger /3. While the peak 
of the ground state is very narrow, the peaks for 
the higher states have larger widths, the reason 
for which is not understood at present. 

3.1. Mass 

The mass of a state is determined from the po- 
sition of the peak in the spectral function. To 
check the reliability of this analysis, the masses of 
the ground and the first excited state obtained in 
this way are compared with those from the double 
exponential fit of propagators for the point and 
smeared sources. At P — 5.90 for p meson, the 
two results for the ground state agree, and those 
for the first excited state are consistent within 
errors, as seen in Fig. ||. 

The first excited state masses for pseudoscalar 
and vector mesons in the chiral limit at each /3 
are extrapolated to the continuum limit in Fig. 
The error for the vector meson at P — 5.90 is 
smaller than that obtained by the exponential 
fit (open square), and our result is consistent with 
that of ref. |S](opcn triangle). In the contin- 
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Figure 2. Comparing the results with these from 
the exponential fit. 
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Figure 3. Continuum extrapolation for the first 
excited state masses for both mesons. 

uum limit we obtain rriT^^ = 660(590) MeV and 
= 1540(570) MeV, which are consistent with 
the experimental values, though the errors are 
large. 

3.2. Decay constant 

The decay constant is related to the area 
of the peak in the spectral function i.e., 
fl cx 2mqJ^^^^,duj fpsiuj)/ml for pseudoscalar me- 



son 



and /2 cx / ^dt^ fv{uj)/ml for vector 



meson. 



As shown in Fig. |J, the ground state decay con- 
stants for both mesons in the chiral limit are con- 
sistent with the results obtained by the exponen- 
tial fit at each /3 and in the continuum limit. 

The decay constants are determined also for the 
first excited state. For the excited tt meson the 
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Figure 6. Divergent continuum extrapolation for 
the unphysical state mass for both mesons. 
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Figure 4. Decay constant for the ground state for 
both mesons. 
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Figure 5. Decay constant for the first excited 
state of p meson. 

decay constant vanishes in the chiral hmit, since 
/tti (X mq/m^.^ by definition. On the other hand, 
in Fig. 1^, the decay constant for p meson is finite 
in this hmit, and we obtain fp-^ — 0.085(36) in 
the continuum hmit. 

3.3. Unphysical state 

We observe in Fig. |l| that there is a peak at 
uja « 2 almost independent of f3. A similar peak 
is also present in the pseudoscalar channel. Ap- 
parently the masses of these states in physical 
units diverge in the continuum limit, as seen in 
Fig. ^ Here we interpret these states as bound 
states of two fermion doublers jsj: In free theory 
the mass of a two doublers system, each with a 
spatial momentum vra, is 2 x log(3) « 2.2, which is 
interestingly close to Lua « 2. The difference may 



be ascribed to its binding energy. We note that 
the bound state of a physical quark and a doubler 
can not appear at zero spatial momentum. 

4. Conclusion 

We have applied the maximum entropy method 
to our precision quenched lattice QCD data 
to extract the spectral functions for pseudoscalar 
and vector mesons. The masses and the decay 
constants extracted from the spectral functions 
for the ground state are consistent with the ones 
determined by the exponential fit, while those for 
the first excited state are new. In future it will 
be interesting to apply the MEM analysis to full 
QCD, QCD at finite temperature and a study of 
decay and scattering. 
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